Gaussian random waves in elastic media 
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Similar to the Berry conjecture of quantum chaos we consider elastic analogue which incorporates 
longitudinal and transverse elastic displacements with corresponding wave vectors. Based on that we 
derive the correlation functions for amplitudes and intensities of elastic displacements. Comparison 
to numerics in a quarter Bunimovich stadium demonstrates excellent agreement. 
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I. INTRODUCTION 

Attracting interest in the field of wave chaos pj , elas- 
tomechanical systems are being studied analytically, nu- 
merically, and experimentally. Weaver first measured 
the few hundred lower eigen frequencies of an aluminum 
block and worked out the spectral statistics @ . Spectral 
statistics coinciding with random matrix theory were ob- 
served in experiments for monocrystalline quartz blocks 
shaped as three-dimensional Sinai billiards [3[, as well 
as, in experimental and numerical studies of flexural 
modes 0,[|| and in-plane modes [f|0] for stadium-shaped 
plates. Statistical properties of eigen functions describing 
standing waves in elastic billiards were first reported by 
Schaadt et al. [§]. The authors measured the displace- 
ment field of several eigen modes of a thin plate shaped as 
a Sinai stadium. Due to a good preservation of up-down 
symmetry in the case of thin plates they dealt with two 
types of modes. The flexural modes with displacement 
perpendicular to the plane of the plate are well described 
by the scalar biharmonic Kirchoff-Love equation [9, 10 . 
In this case a good agreement with theoretical prediction 
for both intensity statistical distribution and intensity 
correlation function was found. However in the case of 
in-plane displacements described by the vectorial Navier- 
Cauchy equation [H, [l(| an agreement between the inten- 
sity correlator experimental data and the theory was not 
achieved Q. 

The aim of present letter is to present an analogue 
of the Berry conjecture for elastic vibrating solids and 
derive the amplitude and intensity correlators with cor- 
responding comparison to numerics. Quite recently, 
Acolzin and Weaver suggested a method to calculate the 
intensity correlator of vibrating elastic solids [ll[ . Based 
on the Green's function averaging technique they suc- 
ceeded to derive the intensity correlator of flexural modes 
generalized due the finite thickness of a plate. Although 
the method might be used for the in-plane modes in elas- 
tic chaotic billiards, we propose here a more simple and 
physically transparent approach based on random super- 
position of traveling plane waves (Gaussian random wave 
(GRW) or the Berry function (T^j). We show that the 
approach allows to derive all kinds of the correlators of 
RGW not only in infinite elastic media but also to take 



into account the double ray splitting at the boundary of 
a plate that plays a significant role in the elastomechan- 
ical chaotic motion [13, LLJ] • We restrict ourselves to the 
two-dimensional case because of the current experiments 
available. Note however that the method can be easily 
generalized for the three-dimensional case. 



II. ANALOGUE OF THE BERRY 
CONJECTURE IN ELASTIC MEDIA 

Shapiro and Goelman [l5[ first presented statistics of 
the eigenfunctions in chaotic quantum billiard although 
their numerical histogram was not compared with the 
Gaussian distribution. This was done by McDonnell and 
Kaufmann [16] who concluded that the majority (~ 90%) 
of the eigenfunctions of the Bunimovich billiard are a 
Gaussian random field. Later it was confirmed by nu- 
merous numerical and experimental studies. The simple 
way to construct RGF is random su per position of par- 
ticular solutions of Eq. (JJ) [l3, G3, [El] with sufficient 
number TV. Thus we come to the Berry conjecture in the 
form [H, [1| 



V'b(x) = 




where the phases n are random distributed uniformly 
in range [0, 2ir) and all the amplitudes are taken to 
be equal (one could assume random independent am- 
plitudes, without any change in the results). The wave 
vectors k„ are uniformly distributed on a d-dimcnsional 
sphere of radius k. It follows now from the central limit 
theorem that both Reybs and ImipB are independent 
Gaussian variables. In a closed billiard the Berry func- 
tion is viewed as a sum of many standing waves, that is 
simply real or imaginary part of function {!]) . 

In our case one has to construct a RGW-function de- 
scribing acoustic in-plane modes. These modes are de- 
scribed by a two-dimensional Navier-Cauchy equation 

/iV 2 u+ (A + /i)V(Vu) +pfl 2 u = (2) 
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where u(x,y) is the displacement field in the plate, A, p, 
are the material dependent Lame coefficients, and p is the 
density. Introducing elastic potentials and A with the 
help of the Helmholtz [2fl| decomposition the displace- 
ment field u could be written, 



u = u/ + u t , u/ = V^, u t = V x A 



(3) 



Eq. © reduces to two Helmholtz equations for the elastic 
potentials 



- v 2 v - fcfV, 

-V 2 A 



fc t 2 A. 



(4) 



Here ki — lj/ci, k t = u>/ct are the wave numbers for 
the longitudinal and transverse waves, respectively and 
uj 2 = pfP/E, where E is Young's modulus. In the two- 
dimensional case potential A has only one none-zero com- 
ponent A z and the dimensionless longitudinal and trans- 
verse sound velocities c;,t are given by 
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where a is Poisson's ratio [10|, 12011 . E and a are func- 
tions of the Lame coefficients [l(j, Our conjecture 
is that both elastic potential be statically independent 
Berry-like functions (JTJ) . We write the potentials in the 
following form 

rr N 

= ^-\/— Vexp[i(k ; „x + 6» ; „)], 

n— 1 

a rr N 

Az(x) = — y — ^exp[i(k tn x + tn )], 



(6) 



n=l 



where 9i n , 9 tn are statistically independent random 
phases. The wave vectors k;„ and k tn are uniformly dis- 
tributed on circles of radii ki and kt respectively. Ac- 
cording to ([3]) the components u, v of the vectorial dis- 
placement field u could be now written 



u(x) 
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(7) 



where </>/„, </>t„ are the angles between k/„, k t „ and the 
x-axis respectively. The prefactors ai = y^y, a t = -y^l — 7 
are chosen from the normalization condition (u^u) = 1. 
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FIG. 1: Correlation function (I10p for 7 = (green solid line), 
7 = 0.5 (red dashed line), and 7=1 (blue dot-dashed line); 
a = 0.345 (aluminum). 



and (...) means average over the random phase ensem- 
bles. Parameter 7 ranges from (pure transverse waves) 
to 1 (pure longitudinal waves). By similar way one can 
construct the elastomechanical GRW for a closed system. 
One can see that the Berry analogue of chaotic displace- 
ments ([7]) is not a sum of two independent GRWs (or 
two independent Berry functions) aitpi + atipt as it was 
conjectured by Schaadt et al. Q with arbitrary coeffi- 
cients a/, a t . In fact, each component u and v in Eq ([7]) 
is related to the Berry functions ([6]) via space derivatives 
in accordance to relations ([3]) . 



III. CORRELATION FUNCTIONS 

First, we calculate the amplitude correlation functions 
in chaotic elastic plate for in-plane GRW ([7]). For quan- 
tum mechanical GRW fl} the two-dimensional correla- 
tion function 

(^b(x + s)Vb(x)) = Ms) 

was found firstly by Berry [T2I ]. Straightforward pro- 
cedure of averaging over ensembles of random phases 
&lm @tn and next, over angles of k- vectors gives 

(it(x + s)u(x)) = ^(cos 2 a f(fas) + sin 2 a g(kis)) + 

-^2(sin 2 a f(k t s) + cos 2 a g(k t s)), 
(u(x + s)«(x)) = 77 (sin 2 a f(kis) + cos 2 a g(kis))+ 

1 ~J I (cos 2 a f(k t s) +sin 2 a g(k t s)), 
(u(x + s)v{x)) = sin 2a(i^ J 2 (k t s) - %J 2 (k lS )), (8) 

where 

f(s) = J (s) - Ms),g(a) = Ms) + Ma). (9) 

It is important to note that the correlation functions © 
were obtained for given direction of the vector s where 
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a is the angle included between vector s and the x-axis. 
However, averaged over all directions of s, the first two 
correlators simplify 



C(s) = (it(x + s)w(x)) = («(x + s)u(x)) 
= Jo(ks) + %J {k t s) 



(10) 



while the third vanishes (u(x + s)w(x)) = 0. One can 
see that in the averaged case the amplitude correlation 
function is defined by two scales because of two different 
sound velocities q, Cj; that is obvious. The correlation 
function C(s) is shown in Fig. [TJ 

Next, we calculate the intensity correlation functions 
P(s) = (7x + s)/(x)) where the intensity I = |u| 2 pro- 
portional to the elastic energy of the in-plane oscilla- 
tions. In quantum mechanics this value is analogous to 
the probability density, the correlation function of which 
was calculated by Prigodin et al. [H|. For the in-plane 
chaotic GRW of the form aitpi + atipt Schaadt et al. [1] 
derived the intensity correlation function as 



P(s) = 1 + 2[af J {k lS ) + a 2 t Jo(hs)f 



(11) 



Our calculations similar to those as for the amplitude 
correlation functions ijHJl give the different result 



P(s) = 1 + ±[("fJo(kis) + (1 - 7 )-/o(M] 2 



+±[{lHk l s)-{l- 1 )J 2 {k t s)f 



(12) 



where r\ = 1 for real GRW and ij — 2 for complex one. 
Although the first term in (fT2|) corresponds to (111|) there 
is a different term consisted of the Bessel functions J2. 
The mathematical origin of deviation is that formula ([7]) 
contains the contributions of the components of the wave 
vectors k; and kt via space derivatives. 



IV. WAVE CONVERSION AT BOUNDARY 

Waves propagate freely inside the billiard, that is, the 
longitudinal and transverse components are decoupled. 
Wave conversion occurs at the boundary according to 
Snell's law 



a sin(0 t ) = c t sin(0;), 



(13) 



The reflection amplitudes for each event of the reflection 
can be easily found following the procedure described in 
[Icj |. At first we consider more easy case of the Dirichlct 
BC (the boundary is fixed). Approximating the bound- 
ary as the straight lines for the wavelengths much less 
than the radius of curvature we have for the reflection 
amplitudes 
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FIG. 2: The intensity I = |u| 2 of the eigen state at frequency 
uj — 28.4 in the quarter of the Bunimovich billiard with fixed 
boundary, a — 0.345. 
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(14) 



Next, we assume that all directions of waves are statis- 
tically equivalent. Then we have for the energy density 
of reflected wave 



Pout = l(T u + T lt ) + (1 - 7)(T« + Ta), (15) 



where 



= 1 
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- / tLdffi, i = l,t. 



Substituting into here l|14[l one can obtain after elemen- 
tary calculations 



// 1 — 777 ^ 1 - Tit — h- 
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T tt = l-f arcsin| + (f 
T ti = |arcsina_fa) 2 7 2 . 
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(16) 



We do not present here integrals Ji, J2, since after sub- 
stitution of (|16p into (|T5|) they cancel each other. The 
equality pi n = 1 = p out gives a very simple evaluation 



(17) 



The next remarkable result is that although the reflection 
amplitudes for the free BC [l(| have the form different 
from (fT4|) (see, for example, formulas in (Ioj|) , the eval- 
uation of 7 by the same procedure gives the same form 
as for the fixed BC. Therefore, we can conclude that the 
result does not depend on either the free BC or the fixed 
BC is applied. Using ([5]) formula (fl"T|) could be written 
in a more simple form 



(18) 
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FIG. 3: Cross marks show numerical results for ratio 7 aver- 
aged over 200 eigenfunctions of the quarter Bunimovich plate. 
The solid line is plotted by formula {TSJ. 



V. NUMERICAL RESULTS AND CONCLUSION 

For numerical tests we took the quarter of the Buni- 
movich billiard and calculated the eigenstates of the 
Navier-Cauchy equation ^ with the fixed BC: u = 
0, v = at the boundary of the billiard by the finite- 
difference method. An example of the eigenstate in the 
form of intensity / = |u| 2 is presented in Fig. O First 
of all we verified formula (fT5|) . For each value of Pois- 
son's ratio a in the range [0, 0.5] with the step 0.05, 200 
eigenfunctions of the billiard were found to calculate av- 



eraged 7. The resulted dependence of 7 on a is shown in 
Fig. [31 that demonstrates a good agreement with formula 
([5]). Thereby we can evaluate 7 for specific a = 0.345 




FIG. 4: The intensity correlation function (|12|l compared to 
the numerics for the same situation as in Fig. [2] 

which correspond to aluminum plate and plot the corre- 
lation functions. The intensity correlation function (|12p 
is shown in Fig. 0] compared to the numerics calculated 
for the eigen function presented in Fig. [2] One can see a 
good coincidence of the theory with the numerical results 
that demonstrates correctness of the GRW-approach ([7]) 
to chaos in elastic billiards. 
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